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Abstract 

We introduce a simple measure of "classicality" of pure and mixed quantum states 
as a maximum value of the Hilbert-Schmidt "scalar products" between the renor- 
malized statistical operators of the state concerned and all displaced thermal states. 
Choosing Fock states as the reference set, we introduce the measure of "anticlassical- 
ity". Both measures are illustrated for the Fock, coherent phase, and generic mixed 
Gaussian states. Gaussian states are shown to be the closest to thermal states pos- 
sessing the same degree of quantum purity. On the contrary, Fock states appear to be 
more close to mixed thermal states than to pure coherent states. 
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1 Introduction 



Since Glauber's paper [1], quantum states for which the so-called P-distribution [2] is non- 
positive or more singular than delta function are called "nonclassical" . For the past decades, 
many authors proposed different quantitative measures of "nonclassicality" . It seems impos- 
sible to reduce all variety of quantum states in infinite-dimensional Hilbert spaces to some 
unique parameter. Therefore, different existing approaches should be considered sooner as 
complementary rather than competitive. 

Historically, the first approach was based on the analysis of deviations from the Poissonian 
photon statistics inherent to the Klauder-Glauber-Sudarshan [1-3] coherent states, which 
are known to be the only "classical" pure states [4-6]. The examples are Mandel's "Q- 
parameter" [7] and its various generahzations [8-10]. The second direction is to evaluate 
the volume of that part of phase space where some quasiprobability distribution assumes 
negative values [11-19] (see also [20,21]). There are also other approaches [22-24]. 

Here we follow the direction opened by Hillery [25]. It consists in evaluating some kinds 
of distances in the Hilbert space between the state concerned and a family of states which 
are assumed to be "classical" . There exist many different definitions of the distance, besides 
the "trace distance" used in [25]. In particular, distances along geodesies on curved manifold 
supplied with Riemannian metrics for coherent, squeezed, displaced, and other states were 
studied in [26-28]. Other examples are the Monge distance [29], "polarized" and "classical- 
like" distances [30] . The most simple from the point of view of calculations is the Hilbert- 
Schmidt distance, used in [30-36], 



Here p is the statistical operator of the quantum state concerned, and pc is related to the 
reference "classical" state. 

For pure quantum states, p = the reference states are usually identified with 

coherent states Pc = \a){a\, as soon as the latter are assumed to be "the most classical" 
ones [4-6]. Then the calculation of distance is reduced to calculating the scalar product 
{ck\i/j). This scalar product is the principal ingredient of many other distances, such as the 
Fubiny-Study distance [37] dpg = 1 — |('(/'i|'(/'2)P (we use here slightly different notation), or 
Wootters' distance [38,39] dw = cos~^(| ('?/^2|V'i) I) (see [30] for more references). 

However, coherent states (whose P-distributions are delta-functions) represent only a 
small subset "at the border" of the set of all "classical" states. Therefore, it seems reasonable 
to enlarge the family of reference states pc, incorporating [35] all displaced thermal states, in 
complete correspondence with original paper [1]. 

Unfortunately, even the Hilbert-Schmidt distance becomes complicated in such a case, 
because the presence of three terms given in the second line of Eq. (1) does not permit to 
find the minimum of this expression analytically (except for the most simple cases), when 
one deals with mixed states. A possibility to simplify calculations by adjusting the "purity" 
Tt{pI) of the reference states to the purity Tr(p^) of the state concerned was discussed in [35]. 

The aim of our article is to consider another possibility, which permits us to find analytical 
expressions for most known families of nonclassical states. Namely, in the next section we 
introduce a new "classicality measure", which is proportional to {but not identical with) 



dnsiP, Pc) 



Tr(p - Pc)' 

Tr(p2) + Tr(p2)-2Tr(pp,). 



(1) 
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the last term in (1), Tr(ppc)- We analyze this measure for the Fock and generic Gaussian 
mixed states. In Sec. 3 we consider the "anticlassicahty measure" , based on using "the most 
nonclassical" Fock states as the reference ones. A brief discussion of results and perspectives 
is given in the concluding Sec. 4 

2 Classicality measure 

As a matter of fact, practically all information about the closeness of the states p and 
Pc is contained in the term Tr(ppc), whereas the presence of the term Tr(pg) in (1) only 
complicates the search of minimum of df^^, adding no significant information. Therefore, it 
seems reasonable to consider only the term Tr(ppc), searching for its maximum with respect 
to the enlarged family of "classical" states, consisting of all displaced thermal states [40-42] 

Pc^D{a)pthD-\a). (2) 

Here D{ol) = exp (^aa^ — a*aj is the well-known displacement operator [1-3] written in terms 
of the bosonic ahhihilation and creation operators, [a, a'] = 1, and pth is the thermal state, 

oo 

Pth^{l-V)j:v'\k){k\, 0<r)<l. (3) 

k=0 

If we confine ourselves to the special case of ?7 = 0, when pc = \a){a\, and consider only 
pure quantum states p = \i!j){ib\. then the quantity max^ |(q;|'?/')P seems a good measure of 
"classicality", because it equals 1 for coherent states and is less than 1 for all other pure 
states. However, dealing with the generic case of mixed states, one meets certain difficulties: 
on the one hand, the measure of classicality for any thermal state must be the same as for 
coherent states, i.e., it must be equal to 1 according to the normalization chosen, but on the 
other hand, Tr(p^) < 1. Therefore, one has to find a suitable generalization of the scalar 
product between pure states to the case of mixed states, which would result in the unit value 
for identical states. 

One possible solution, satisfying many additional requirements, was given by Uhlmann 
[43], who showed that the generalization of the quantity |('0i|'02)| can be taken in the form 

Tt\/ Px"^ P2P\'^ . However, the calculation of this trace is rather complicated problem, al- 
though certain progress has been achieved recently for Gaussian states (in the case of infinite- 
dimensional Hilbert space) [44-47]. 

We prefere to follow a more straightforward way: simply to replace the operators p and 
Pc in Tr(ppc) by the renormalized operators p' and p'^ , where 

p'^p/^T^). (4) 
Therefore, we define the "classicality measure" as 

C = maxIV(p'p'J. (5) 

Pc 

For pure states, p = we have 

C = m^ax(^|p;iV'). (6) 



3 



In modem literature, the "generalized scalar products", such as Tryp^ P2P1 or Tr(/3i/)2), 
are frequently called [48,49] "fidelities" ( "Bures-Uhlmann" or "Hilbert-Schmidt" , respec- 
tively). In particular, a possibility of using the "modified fidelity" Tr(p'^p2) was mentioned 
in [50]. Note that if only one of the states is pure, then the "modified fidelity" does not 
coincide with {ipi\p2\'4'i) but differs from it by the factor [Tr(p2)]~^/^. Sometimes [48,49] this 
is considered as a flaw. But in our problem the factor discussed has real physical meaning, 
because it permits us to distinguish mixed states from their pure "partners" possessing the 
same probability distribution functions of quanta (see examples below). Another difference 
between the cited papers and our approach is that we consider "fidelity" not to a fixed state, 
but to the whole family of "classical" states. 

We shall use the notation /(?7,q;) for the function Tr(/9'/3^), whose maximum gives the 
value of classicality. Sometimes, it can be difficult to find a maximum with respect to both 
variables. Then one can calculate partial maxima with respect to one variable, fixing zero 
value for another. We shall denote such "reduced classicalities" as Ca or C^, where the 
subscript indicates the parameter, with respect to which the maximum was calculated. 



2.1 Classicality of the Fock states 

For the Fock state the quantity Tr(ppc) is reduced to the probability {n\pc\n) of dis- 
covering n quanta in the displaced thermal state. This probability was calculated by many 
authors, e.g., in [40-42,51]. Using the results of [51] and taking into account the expression 
for purity of the thermal state (it does not depend on the shift parameter a), 

we arrive at the function 



fn{V:Oi) = y'l - ?]2r^"exp (^-|ap^l - 772 

xL„(-H^(l- 77)777), (8) 

where Ln{z) is the Laguerre polynomial. For 7/ = we have the well-known Poissonian 
distribution of the coherent state, 

/n(0, a) = exp (-|q;|^) {af^'/nl, 

which has maximum at [ap — n. Therefore, confining ourselves to pure coherent reference 
states, we obtain for the reduced classicality of the Fock state the formula (cf. [35]) — 
e~^n'^/n\, which results in Cj^^ ~ (2™)"^/^ for n 1. 

On the other hand, the function fniv^ 0) has maximum at 77^ = ^Jn/{n-\- 1), so that 



e^^^ = ^n^/{l+nY+K (9) 

One can verify that Cj") > C'^^'^ for any n > 1. In particular, = 1/2 and C^^^ = 1/e, 
whereas for n » 1 we have (en)"^/^. Moreover, the function 7^(77, a) has negative 



4 



derivative with respect to |ap at the point lap = for any rjn- Therefore, the maxima 
of fn{ri,a) at the points (r^n, 0) are global, whereas the maxima at the points (0,n) are 
only local: see Fig. 1. This means that the classicality of the Fock state is given by (9), 
and that the pure state \n) is more close not to the pure (displaced) coherent state with 
Ucofi = lap = n, but to the mixed (undisplaced) thermal state with 

rith = 7]/{l-rj) = y/n (^^/n + Vn + l) 

(i.e., nth ~ 2n for n ^ 1) and 

Tr(p^) = (V^+V^)"'. 



2.2 Classicality of Gaussian states 

Generic Gaussian states are characterized [51] by two variances of the quadrature compo- 
nents, (Tq and (Tp, their covariance, apg, and two displacement parameters in the phase plane 
qp. It is clear, however, that for the fixed (co)variances, the quantity Tr(ppc) is maximal for 
coinciding displacement parameters of p and pc, which means that it is sufficient to consider 
the states with zero displacements, i.e., to find the maximum of the function of single vari- 
able F{ri) = f{r],0). Since the statistical operator of the thermal state is diagonal in the 
Fock basis, the function F{ri) can be easily expressed in terms of the generating function 

oo 
n=0 

for the photon distribution p„ = {n\p\n): 

F(n)^p-'/'^fl^^G{rj), /.^TV(p^). (10) 

The function G{z) for generic Gaussian states was calculated by means of different methods, 
e.g., in [41,51,52]. Using the form given in [51], we arrive at the function 



F{rj) = 2 



a — 2br] + erf 



1/2 



(11) 



where 



1 + At' + 2Tp\ h^l-n\ 

-1 



4 (cTpCr, 



<1 ^pq 



c = 1 + /x^ - 27//^ 

T — (Tp-\- Uq. 



One can easily verify that the maximum of function (11) is attained for 77^. = (1 — /x)/(l + /i), 
which means that the thermal state which is most close to the given Gaussian state has the 
same purity: pLc — A*- The classicality of the Gaussian states is given by a simple formula 



l + p,T 



(12) 
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Note that jjT > 1 for Gaussian states, and the equahty holds for thermal states. The 
parameter T is related to the mean number of quanta n in the unshifted Gaussian state 
as T = 1 + 2n. Therefore, the classicality of pure unshifted Gaussian states (which are 
nothing but squeezed vacuum states) equals C^*^"^ = (n + If < 1, then the minimal 

mean number of quanta in the Gaussian state equals n„tm = (1 ~ /^)/ (2/^)- However, mixed 
Gaussian states remain unsqueezed (that means, no one quadrature component can have 
the variance less than the vacuum state value 1/2 in dimensionless units), ii n < Tie — 
(1 — /x^)/(2//^). Such states possess positive P-distributions [40], and it was proposed in [47] 
to extend the family of "classical" states, including all mixed unsqueezed Gaussian states. 
We suppose to study such an "extended classicality" in another paper. 

There exist pure states possessing the same photon distribution function as thermal ones; 
usually they are called coherent phase states [53, 54] : 



k> = \/i-l^l'E^"l^)- (13) 

n=0 

Although the state (13) is most close (at least with respect to the reduced C^-measure) to 
the thermal state with r) — jsj^, its r^-classicality is less than 1 (precisely because we use 
renormalizcd operators in the definition (5), the function F{r]) in Eq. (10) contains the 
factor /i"^''^, but this factor equals unity now, as soon as p'^ = p^): 



V 



\ 



1 - 


e 


2 


1 + 


£ 


2 



[1 + 2n) 



in accordance with our feeling that pure states, in general, are less classical than their 
mixed partners. For highly excited states with the same mean value n ^ 1 we obtain the 
asymptotical relations 



3 Anticlassicality measures 

As soon as the Fock states are usually considered as "the most quantum" states, we can use 
them as reference states for the definition of the degree of "anticlassicality" : 

^ = max(n|p|n). (14) 

Perhaps, it could be worth including the displaced Fock states, \h) = D{a)\n), into the set of 
reference states, but here we shall not do it. Actually, we can define two measures: one given 
by (14), where integer n runs over all integers, including n = 0, and another (denoted as ^i), 
where the vacuum state |0) (which is distinguished from all other states) is excluded from 
the set of reference states. To see the difference, let us consider an example of coherent states 
\a). One can easily find that the dependence of A on the mean photon number n — \a\'^ has 
different analytical forms in the intervals k <n < k + 1, k = 0,1, . . . : 

^H(ri) =exp(-ri)nV^!, k<n<k + l. (15) 
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The function given by (15) is continuous, but it has jumps of derivatives at the points rik = k: 
the right derivatives at these points are equal to zero, while the left derivatives have finite 
negative values. However, A^°'\n) goes to 1 as ri — > 0, and this may seem strange. 

The function Ai'\n) coincides with A''°'\n) ior n > 1. But for n < 1 its behaviour is 
quite different: 

A^{'\n) = n exp (— n) , < n < 1. 

The dependence A^"\n) is shown in Fig. 2. Asymptotically, for n ^ 1, we have (using the 
Stirling formula) A^i"^ ^ {27in)^^^^. 

As another example we consider the squeezed vacuum state. It is known that the photon 
distribution function p„ of this state is different from zero only for even values n — 2m, and 
it can be expressed in terms of the mean number of photons as [51] 

-^1 + ^^ (2-m!)Hl+J ■ ^^^^ 
The sequence p2m decreases monotonously. Therefore we obtain a monotonous dependence 

whereas the function _ 

^""^ 2(l + n)3/2 

has maximum ^[^11 = l/(3\/3) for n = 2 (and the most close Fock state is |2)). 
For the coherent phase state (13) wc obtain 

A^'\n) = {l + n)-\ (n) =71(1 +71) 

with A^i^ax — 1/4 at ri = 1. According to definition (14), the same results must hold for 
thermal states, since they have the same photon distribution function {n\f)\n). However, 
we can distinguish pure and mixed states having the same diagonal elements in the Fock 
basis, playing again with the purity parameter ii = Tr(p^). But in contrast to the case of 
classicality parameter, now we should, sooner, multiply the quantity {n\p\n) by /x (or some 
power of II, if one wishes), rather than divide it by /i (or s/JI), because the intuition tells us 
that mixed states are more anticlassical than pure ones. Adopting the definitions 

A — /J. m.a,'x{n\p\n), Ai — p m.a'x{n\p\n) , (17) 

we obtain for the thermal states the expressions 

A^^^^n) = [(1+7I)(1 + 27I)]-' 



Ann) 



;i + n)2(l + 2n)' 

These three cases are illustrated in Fig. 3. The asymptotical dependences at 7i ^ 1 are as 
follows: 



2v^' ' 2n' ' 2rf 

For generic (mixed and shifted) Gaussian states, the dependences Ai{n) and Ai{n) may 
be more complicated, especially in the regime of strong irregular oscillations of the photon 
distribution function p„ [55,56]. We reserve this case for another study. 
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4 Conclusion 



We have introduced new measures of "classicality" and "anticlassicality" of quantum states, 
which characterize the "fidehties" of the chosen state to the whole famihes of states con- 
sidered as "classical" or "maximally quantum". Due to simplicity, these measures can be 
calculated analytically for many important sets of quantum states, thus providing additional 
information on their properties. In particular, we have shown that the most close to arbi- 
trary Gaussian states are thermal states with the same degree of quantum mixing (purity). 
At the same time, there are (mixed) thermal states which are more close to the Fock states 
than (pure) coherent states. As soon as the new measures distinguish coherent quantum 
superpositions ("cat states") from quantum mixtures, they could be used in studies of de- 
coherence processes. Also, it would be interesting to establish relations between the new 
and existing measures of classicality, such as Mandel's parameter, "nonclassical depth" [13], 
and so on, or to try to "order" different families of quantum states in accordance with their 
degrees of (anti) classicality. We hope to report on the results of studies in these directions 
somewhere. 
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Figure 1: The renormalized Hilbert-Schmidt fidelity between the Fock state |3) and displaced 
thermal states. 



Figure 2: The anticlassicahty measure Ai'^ of coherent states versus the mean number of 
quanta n. 



Figure 3: The anticlassicahty measures Ai{ri) for the squeezed vacuum and phase coherent 
states, compared with the modified measure Ai{n) for the thermal states. 
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